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ABSTRACT 

A generalized one dimensional computer code 
for analyzing the flow and heat transfer in the tur- 
bomachinery cooling passages has been developed. 
This code is capable of handling rotating cooling 
passages with turbulators (inline, staggered, or in- 
clined), 180 degree turns (both sharp and gradual 
turns, with and without tur oulators) , pin fins, finned 
passages, by-pass flows, tip cap impingement flows, 
and flow branching. The code is an extension of a 
one-dimensional code developed by P. Meitner. In 
the subject code, correlations for both heat transfer 
coefficient and pressure loss computations have been 
developed to model each of the above mentioned type 
of coolant passages. The code has the capability of 
independently computing the friction factor and heat 
transfer coefficient on each side of a rectangular pas- 
sage. Either the mass flow at inlet to the channel or 
the exit plane pressure can be specified. For a spec- 
ified inlet total temperature, inlet total pressure and 
exit static pressure, the code computes the flow rates 
through the main branch and the sub branches (if 
any), flow through tip cap for impingement cooling 
(if any), in addition to computing the coolant pres- 
sure, temperature and heat transfer coefficient distri- 
bution in each coolant flow branch. Predictions from 
the subject code for both non-rotating and rotating 
passages agree well with experimental data. The code 
Wets used to analyze the coolant passage of a research 
cooled radial rotor. 

INTRODUCTION 

Increasing the turbine inlet temperature is a 
very lucrative approach to improving the thermal ef- 
ficiency and specific thrust of gas turbine engines. 
However, this practice is limited by the availability 
of materials which can withstand these high temper- 
atures. Current materials cannot be exposed to tem- 
peratures above 1250 K without a dramatic decrease 
in life [1]. Blades must be cooled if significantly 


higher turbine inlet temperatures are desired. As the 
cooling scheme becomes more complex, a better un- 
derstanding of the internal pressure and heat transfer 
coefficient distributions is necessary to be able to pre- 
dict the hot section metal temperatures accurately. 
Typical cooling concepts of modern turbine blading 
are shown in figure 1. Figure 1(a) shows the cooling 
scheme for an axial rotor, while figure 1(b) shows a 
hypothetical cooled radial rotor. The shown cooling 
passages may have heat transfer enhancement devices 
in the form of turbulators, pin fins, finned passages, 
leading edge tip cap impingement flows, etc. The 
main coolant flow may also be split into sub branches 
as shown in figure 1(a) and 1(c). It should be noted 
that to obtain a 2-D distribution of heat transfer co- 
efficient in the plane of the blade surface will really 
require a 3-D analysis. Because of the complex nature 
of the flow field in cooling passages, it is extremely dif- 
ficult to develop a generalized 3-D flow analysis code 
capable of modeling different cooling passage config- 
urations and obtaining detailed distribution of heat 
transfer coefficients and pressure drops. Such infor- 
mation is required in the detailed thermal and struc- 
tural analysis of the rotor and blading [2]. Hence it 
was felt that development of a generalized 1-D code 
which can model various types of cooling passages in- 
cluding flow branching should be useful in the design 
and analysis of cooled turbomachinery. Based on the 
literature review, the authors conclude that such a 
generalized 1-D code is not available in the open lit- 
erature. 

The code reported in this paper is an extension 
of the code developed by Meitner [3], and is appli- 
cable to any combination of cooling passage config- 
urations mentioned above. Specifically, the reported 
code differs from the code of reference [3j in that it 
can calculate and balance flows in sub-branches, and 
can thus compute flows in passages shown in figure 



1(a) and 1(c) whereas the code of reference [3] will 
not be capable of handling these passages. Also, 
the code reported in this paper can make separate 
calculations for the leading(suction) surface, trail- 
ing (pressure) surface, and side walls of a rotating 
channel of rectangular cross section. This is very 
important for capturing the experimentally observed 
behavior in rotating environment. In addition, a 
more detailed modeling of each type of cooling pas- 
sage above has been incorporated in the present code 
compared to the code of Meitner. The code can pre- 
dict the mass Dow rates, internal pressure and tem- 
perature distribution, heat transfer coefficient and 
heat Dux distribution at the walls along the passage. 
In a limited sense, it can predict the heat transfer 
coefficient distribution along each side of rectangu- 
lar passage. It integrates the 1-D compressible mo- 
mentum and energy equations along a defined flow 
path and will account for area changes, mass ad- 
dition/subtraction, Coriolis acceleration, centrifugal 
force effects, bouyancy force effects, friction and heat 
addition. Flow can be bled off for tip cap impinge- 
ment cooling and a flow bypass* can be specified in 
which the coolant flow is taken ofT at one point in 
the flow channel and re-introduced at a point fur- 
ther downstream in the same channel. As mentioned 
above, other types of passages that can be modeled 
include plain passages, rib-roughened passages (stag- 
gered and inline turbulators at various angles of at- 
tack), different pin fin arrangements (with or without 
lateral flow ejection, as shown in figure 1(a)), finned 
passages, 180 degree turns (both sharp and gradual) 
and flow branching. 

Based on the most recent experimental data 
available, correlations have been obtained to model 
the various types of rotating cooling passages. The 
user has the option of specifying the temperature and 
pressure at inlet to the main branch and either the 
main flow rate or the exit pressures from the sub 
branches. The present version of the code cannot 
accomodate supersonic flow anywhere in the cooling 
passage. The flow can at best be choked at exit from 
one or more of the branches. However, it is felt that 
this is not a serious restriction for modeling turboma- 
chinery cooling passages. 

Comparisons have been made with a wide vari- 
ety of non-rotating and rotating passage experimen- 
tal data representing various types and combinations 
of cooling passages. It is shown that the predictions 
compare well with the experimental data. 



(a): Axial Turbine [9] 



(b): Radial Turbine 



(c): Radial tubine with Flow branching 

Figure 1. Typical Cooling Concepts of Modern Gas 
Turbines 
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FLOW PATH MODELING AND 
COMPUTATION PROCEDURE 

The flow path is defined in cartesian co-ordinates 
z,r and y as shown in figure 1(b). The flow path ge- 
ometry is described at nodes and the type of passage 
between each set of nodes is specified using different 
parameters in the code input. 

For specified entrance conditions and a specified 
(or guessed) flow rate, the flow solution is obtained 
by marching along the defined flow path from node 
to node. The pressure and temperature at each node 
is calculated from the values at the previous node by: 



( 4 ) 


Pn = Pn-l + ( dp/dx) n * dx n (1) 

T n — T n -1 + (dT/dx) n * dx n (2) 

where (dp/dx)n and (dT/dx)n are given by the solu- 
tion of momentum and energy equations as [3]: 


lm/ dx) n = (m„ - m n -i)/(l n ) 

(5) 

(dA/dx) n = (A n - An- l)/(ln) 

(6) 

(dr/dx) n = (r n - r n -i)/(l n ) 

(7) 
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The constants Cl, C2 and CT are defined as fol- 
lows [3]: 


Ci = 


(m/A) 2 R 
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( 8 ) 


_ (m/A) 2 RT Ci 

°2 - 1 -X + -=■ 


pC P 


{m/A^RT 2 


(9) 


Ct = 1 + 


(m/A) 2 R 2 T 

P 2 C p 


(10) 


Iteration takes place in each marching step until 
a convergence is acheived. After marching through all 
the nodes, the calculated branch or sub branch exit 

f pressures are compared with the given exit pressures 
if specified). If necessary, a correction is made to the 
inlet flow rate and the marching procedure repeated 
until the calculated and specified exit pressure values 
agree within a prescribed tolerance. 

For flows where the local coolant passage Mach 
number becomes high (> 0.9), pressure changes in the 
flow direction become highly non-linear. A provision 
has been made in the code to automatically add ad- 
ditional nodes between the original nodes and repeat 
the iteration process mentioned above until the Mach 
number values at each of the original nodes do not 
change between successive increases in the number of 
nodes. 

In the present paper, for the different types of 
cooling passsages mentioned above, both friction fac- 
tor and heat transfer coefficients are modeled with 
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correlations developed from experimental data avail- 
able in the literature or using the correlations avail- 
able in the literature (with modifications). For val- 
idating the correlations developed, computations us- 
ing these correlations have been compared with ex- 
perimental data from an entirely different source than 
those used in the development of correlations. 

The correlations developed are generally of the 
form: 


Nu or f = a(NDl) b (ND2) c (ND3) d ... (11) 

where Nu is the Nusselt number, f is the friction 
factor, ND1, ND2, etc are relevant non-dimensional 
quantities like Reynolds number (Re), Rotational 
Rayleigh number (Ra), Rossby number (Rob Prandtl 
number (Pr), rib height (or roughness height) to hy- 
draulic diameter ratio (e/J9^), the rib pitch-to-height 
ratio (P/e), channel aspect ratio (W/H), pin-fin span- 
wise spacing to diameter ratio (S/d), pin-fin stream- 
wise spacing to diameter ratio (X/a), pin height to 
pin diameter ratio (H/d), etc. The correlations were 
obtained from the experimental data using the SAS 
[4] multiple correlation and regression analysis pro- 
cedure. The correlation equations are not shown in 
the paper due to space limitations. Additional de- 
tails of the correlation equations and the code will be 
presented in a future NASA report. 

In the present version of the 1-D code, correla- 
tions have been included to model any combination 
of the following types of passages: 

a) Plain passages (Smooth and Rough) : Based on 

the magnitude of relative roughness parameter 
(e/Dh), the passages are classified as smooth or rough 
and the corresponding correlation equations are used. 

Friction Factor: For laminar flow in circular, rectan- 
gular or square ducts, based on the data available in 
Kays and Crawford [5] in the form of tables and fig- 
ures, correlations have been obtained to model both 
smooth and rough non-rotating passages and for both 
fully developed and entrance regions of flow. For tur- 
bulent flow and in the transition region in ducts, cor- 
relations listed in Welty, Wicks and Wilson [6] have 
been used for both smooth and rough passages 

Heat Transfer Coefficient: For laminar flow, correla- 
tions are developed based on the data of Kays and 
Crawford [5] for the fully developed as well as the 
entrance regions of flows in both circular and non- 
circular ducts. For fully developed turbulent flows, 
correlations given by Kreith and Bohn [7] and/or 
those in reference [5] have been included. Entrance 
effects are modeled based on Kays and Crawford [5]. 

b) Turbulated Passages : Correlations have been in- 
corporated to model the heat transfer coefficient and 
pressure drop distributions for rectangular or circular 
passages with transverse rib rougheners (both inline 
and staggered arrangements and at any given angle 
of attack). Correlations have been included for both 
the turbulated walls and the smooth side walls based 


on the work of Han, et al. [8,9] for in line rib ar- 
rangement and Taslim and Spring [10] for staggered 
arrangement. Correlations include the effects of non 
dimensional parameters like rib height to hydraulic di- 
ameter ratio ( e/Dh ), rib pitch-to-height ratio (P/e), 
channel aspect ratio (W/H), Reynolds number (Re) 
and rib angle of attack ratio (a/90). Based on Boyle’s 
data [11], entrance effects have also been incorporated 
in the correlations. 

c) Pin Fins : Based on the work of Metzger, et al. 
[12,13,14,15], correlations are available to model pres- 
sure drop and heat transfer coefficients (array aver- 
aged) for different staggered pin fin configurations. 
Nusselt number (Nu) correlations which include the 
effects of streamwise pin spacing are given by Metzger 

[12] . The local row resolved heat transfer coefficient 
correction has been correlated based on the data of 
Metzger [12] as well as the corrections to Nu to in- 
clude effects of channel convergence. Nusselt number 
correlations are also available based on the work of 
VanFossen [16], It should be noted that in VanFos- 
sen correlations, the Nusselt numbers are based on 
hydraulic diameter DV which depends on the open 
volume in the pin fin bank and the total heat trans- 
fer surface area whereas the Metzger correlations are 
based on pin diameter. In a recent review paper on 
staggered array pin fin heat tansfer, Armstrong and 
Winstanley [17] recommend the use of Metzger or 
VanFossen correlations for short pin fins (H/d <= 3) 
and the correlation of Faulkner [18] for long pin fins 
(H/d > 3). Friction factor correlations for both short 
and long pin fins are based on the work of Metzger 

[13] - 

The effect of channel convergence on friction fac- 
tor is incorporated based on the work of Brown [19]. 
For pin fin channels with lateral ejection, the correla- 
tions developed by Lau et al. [20] are used. 

d) 180 Degree Turns : Based on the experimental data 
of Metzger, et al. [21,22,23,24], correlations were de- 
veloped for heat transfer coefficient and pressure loss 
in smooth and turbulated sharp 180 degree turns. For 
gradual 180 degree turns, the pressure loss coefficients 
were based on the equations from reference [25] and 
the correlations for the heat transfer coefficient vari- 
ation was based on the experimental data of Baughn, 
et al. [26]. 

e) Finned Passages : For cooling passage channels 

with fins parallel to the flow direction (usually near 
the trailing edges of blades, figure 1(b)), an effective 
heat transfer coefficient is obtained [3] which can be 
applied to the corresponding plain passage surface 
area without fins. 

f) Flow Branching : Pressure loss due to flow branch- 
ing has been modeled based on reference [25]. 

g) Leading Edge Tip Cap Impingement Cooling : 
Modeling of the tip cap impingement cooling is based 
on the development by Meitner [3] who utilized the 
data of Chupp et al. |27]. It includes the effects of 
rotation on pressure changes across the tip cap im- 
pingement chamber and the temperature increase due 
to impingement in the chamber. The model also com- 
putes the mass flow through the impingement holes. 
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Corrections have been applied to the Meitner model 
based on the work of Epstein et al. [28] which include 
the effects of rotation on heat transfer in a impinge- 
ment cooled turbine blade. Modeling of by-pass flows 
(see figure 1 (b)) is also based on Meitner’s develop- 
ment [3], 

h) Rotational Effects of Plain Passages : Effects on 

both pressure drop and heat transfer coefficient due 
to Coriolis forces, centrifugal forces, and centripetal 
buoyancy have been included for ducts (circular, 
square or rectangular) rotating in orthogonal and/or 
parallel modes (with respect to the axis of rotation) 
and for both laminar and turbulent flows. 

For laminar flow in plain passages with smooth 
or rough walls, and with their axes normal to the axis 
of rotation, the pressure loss correlations are based on 
the work of Mori and Nakayama [29] and heat trans- 
fer correlations are based on the work of Morris [30] 
and for moderate temperature differences on the an- 
alytical work of [29]. For turbulent flows, the heat 
transfer correlations are based on the work of Clifford 
[31] or Morris and Ayhan [32,33] with corrections for 
the leading and trailing surfaces based on Harsagama 
and Morris [34]. The user also has the choice of bas- 
ing it on the theoretical work of Mori, Fukada and 
Nakayama [35]. The pressure loss correlations are ob- 
tained from the work of Ito and Nanbu [36] or from 
reference [35]. 

For laminar flow in plain passages rotating about 
a parallel axis, the pressure drop correlations are 
based on the rotating friction factor given by 

/rot = f.t + fex (12) 

f ex = 0.503 Ro 106 (13) 

where f ex is “the excess friction factor” suggested by 
Morris [36]. The heat transfer correlations are based 
on the work of Morris and Dias [38] and for moder- 
ate temperature differences on tne theoretical work 
of Mori and Nakayama [39]. For turbulent flows, the 
heat transfer correlations are based on the work of 
Morris and Dias [40] and Woods and Morris [41] and 
for moderate temperature differences on the theoret- 
ical work of Nakayama [42], The pressure drop corre- 
lations are based on [43]. For the entrance region, the 
heat transfer correlations are based on work of Morris 
and Woods [43]. 

It has been shown from the experimental data 
of Hajek et al. [44], and Guides [45] that the heat 
transfer characteristics due to rotation are radically 
different for the leading and trailing surfaces and also 
for the centrifugal (radially outward) and centripetal 
(radially inward) passages. These features have been 
captured in the correlations used in the present code. 

i) Rotating Channels with Heat Transfer 

Enhancement Devices: 


transfer coefficient and friction factor (pressure 
drop) distribution due to enhancement devices and 
then computing the effects of rotation as a ra- 
tio of a) Nu(rotating)/Nu(non-rotating) and b) 
f (rotating) /f(non-rotating) based on the experimen- 
tal correlations mentioned above. The validity of the 
correlations used is shown in the following section. 

It should be noted that the user always has the 
option of providing the pressure loss coefficient, fric- 
tion factor, and/or heat transfer coefficient for any 
given interval in the form of data to be curve fitted, 
equations, or a set of given values. 

COMPARISON WITH 
EXPERIMENTAL DATA 

In the following paragraphs, the predictons from 
the code are compared with a wide variety of exper- 
imental data for various types of heat transfer en- 
hancement devices. As mentioned earlier, the predic- 
tions from the code are compared with the data from 
the sources which are different from those used in the 
correlation. Both non-rotating and rotating passage 
data with and without various heat transfer enhance- 
ment devices are compared to give a clear picture of 
the range of application of the code. 

a) Non-Rotating Passage Data : 

For a straight duct with 90 degree turbulators, 
the code predictions are compared with the experi- 
mental data of Taslim and Spring [46,47] in figures 
2 and 3. Figure 2 compares the mean Nusselt num- 
ber variation with Reynolds number (Re) for a fixed 
rib height to channel hydraulic diameter ratio (e/D/j, 
also known as blockage factor), rib pitch to height ra- 
tio (P/e) and channel aspect ratio (W/H). The defini- 
tion of these terms are clarified in the inserts in figures 
2 and 3. Figure 3 compares the Nusselt number en- 
hancement, expressed as (Nu /Nu 8t ) f with a variation 
in the blockage factor e/D h for fixed values of P/e, 
W/H, and Re. Predictions are in excellent agreement 
with the experimental data. 

In figures 4 and 5, the code predictions for a two- 
pass serpentine passage with a 180 degree turn, for 
both smooth and turbulated walls, are compared with 
experimental measurements. The calculated Nusselt 
number variation along the passage is compared with 
the data of Chandra [48] in figure 4 and with the 
data of Boyle [11] in figure 5. Agreement with data 
of [48] is very good for both smooth and turbulated 
passages. The reason for the reduction in the Nusselt 
number in the 180 degree turn region ( X/D ^ = 12 
to 14) for the turbulated case is that there are no 
turbulators in the turn region. In the case of Boyle’s 
data, the 180 degree turn region (X/Dh = 14 to 16) 
has turbulators. 

Predictions agree very well with Boyle’s data [ll] 
at higher values of Re where as at low values of Re (< 
= 20000), the agreement is good except in the region 


The combined effects of rotation and heat 
transfer enhancement devices (like ribs, pin fins 
etc.) are obtained by first computing the heat 
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Figure 2. Comparison of Nussclt number variation 
with Reynolds number with the turbulated passage 
data of Taslim and Spring [46,47] 



Figure 3. Comparison of Nusselt number variation 
with blockage factor (e/Zl/,) with the tubulated pas- 
sage data of Taslim ana Spring [46,47]) 



Non-dimensional distance along flow path, X/D 

h 


Figure 4. Comparison of variation of Nusselt number 
with X/Dh with the two-pass rib roughened passage 
(at 90° with a sharp 180 deg. turn) data of Chandra 
[48], 



Non-dimensional distance along flow path, X/D 

n 


Figure 5. Comparison of variation of Nusselt number 
with X/Dh with the two-pass turbulated passage data 
of Boyle [11]. 

immediately after the 180 degree turn where the com- 
puted Nu drops sharply compared to the experimen- 
tal data. The maximum error is of the order of 20 to 
30% in this region. The reason for this behavior is 
not known at present. 

For passages with pin fins, the code predictions 
are compared with the experimental data of Lau et 
al. [49] and Peng [50] in figures 6 and 7. Figure 6 
compares the predicted Nu variation with row number 
with the data of Lau [49] and figure 7 compares the 
predicted variation of friction factor (f) and Colburn 
factor (j) with Reynolds number with the data of Peng 
[50]. Both the figures show that the predictions are 
in excellent agreement with the experimental data. 
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Figure 8. Schematic of the rotating serpentine pas- 
sage of Hajek [44] showing the four zones compared 
with the code predictions. 


Figure 6. Comparison of variation of Nusselt number 
with row number with the pin fin data of Lau [49]. 



Fig me 7 Comparison of variation of friction factor 
and Coiuuwi lactoi with Reynolds number with the 
pin fin data of Peng [50]. 

b) Rotating Passage Data : Extensive comparison has 
been made with the rotating passage data (both 
smooth and turbulated serpentine passages) of Hajek, 
et al. [44]. A schematic of the serpentine passage of 
Hajek is shown in figure 8. A 9 shown in figure 8, four 
typical zones have been chosen for comparison with 
code predictions to illustrate the validity of the corre- 
lations developed for rotating channels. Zones 1 and 
2 are in the radially outward (centrifugal) section and 
zones 3 and 4 are in the radially inward (centripetal) 
section of the passage. As explained earlier, separate 
correlations are used to predict the effects of rotation 
on the leading surface (suction surface!, trailing sur- 
face (pressure surface) and the two side walls of the 
coolant passage. This enables the code to compute 
the rotational effects more accurately. 


Figure 9(a)-(d) compare the Nusselt number 
variation with Rossby number (Ro) for the above 
mentioned four zones along the passage. As seen in 
figures 9(a) and 9(b), the predictions clearly show the 
experimentally observed trend for the centrifugal sec- 
tions zones 1 and 2, that as the Rossby number in- 
creases, due to increasing rotational speed, Nusselt 
number increases on the trailing (pressure) surface 
and decreases on the leading (suction) surface. This 
trend is observed for both smooth and turbulated pas- 
sages. This is due to the increased secondary flow 
effects at higher values of Ro. Figures 9(c) and 9(d) 
compare Nusselt number variation with Ro for the 
centripetal channel. These figures also show that the 
trend of Nu variation with Ro is in general opposite to 
that of the centrifugal section. Comparison of pres- 
sure variations along the passage with the data of 
Hajek [44] is made in figure 10 for the turbulated pas- 
sages. All the above figures clearly show that the code 
predictions agree well with the data of Hajek [44]. 



(a): Zone 1. 

Figure 9. Comparison of Nusselt number variation 
witli Rossby Number for the smooth and 45 deg. 
turbulated serpentine rotating passage data of Hajek 
(44j. 
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Rossby Number, Ro 


(b): Zone 2. 



(c): Zone 3. 
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Figure 10. Comparison of pressue distribution along 
the passage with the turbulated rotating serpentine 
passage data of llajek [44). 

Comparison of Nusselt number variation with the 
centrifugal passage data of Guide* [45] is made in fig- 
ure 11 and with the data of Iskakov and Tfushin [51] 
figure 12. These figures show that the code predic- 
tions are in good agreement with the experimental 
data in botli cases with a maximum error of about 
20%. The trends shown in these figures are the same 
as the centrifugal passage data of llajek [44]. 



Rossby Number, Ro 


Figure 1 1 .Comparison of Nusselt number variation 
with Rossby Number with the data of Guides [45]. 
at X/D h = 7.4. 


(d): Zone 4. 
Figure 9. Concluded 
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Figure 12. Comparison of Nusselt number with 
the parameter Gr/Re with the data of Iskakov and 
lYushin (51j at X/Du = 7* 

COOLED RADIAL TURBINE ROTOR ANALYSIS 

The coolant flow of a cooled NASA research ra- 
dial turbine rotor, which will be tested at the Lewis 
Research Center, has been analyzed using the code. 
The cooling passage consists of a main passage split- 
ting into three branches and includes pin fins and 
fumed passages. The schematic of the cooling pas- 
sage is shown in figure 1(c). For a given inlet total 
temperature and pressure and exit pressures from the 
branches, the present code was used to predict the dis- 
tribution of pressure along the entire passage and also 
the mass flow through the main and sub branches. 
The results are shown in figure 13 and table 1 respec- 
tively. These predictions will be compared with the 
experimental results as they become available. 


Table 1. Mass flow computations for the coolant pas- 
sage of a NASA research radial turbine 


Branch 

m/m, ref 

Code 

Computation 

Main 

6.73d 

B1 

2.399 

B2 

1.632 

B3 

2.703 


m r «/ = 0.00lA:g/s 



Figure 13. Pressure distribution along the coolant 
passage of a NASA researcli radial turbine. 

CONCLUDING REMARKS 

A generalized one dimensional flow code appli- 
cable to compressible flow through complex rotat- 
ing turbine cooling passages with various heat trans- 
fer enhancement devices lias been developed. Even 
though this is a 1-D code, the effects of rotation on 
heat transfer is taken into account in a quasi 2-D fash- 
ion, i.e. a separate correlation is available for eacli 
side of a rectangular passage. The code has been val- 
idated using a wide variety of both non rotating and 
rotating passage data. The authors feel that the code 
will be useful to the turbine designer in evaluating 
the effects of various cooling passage designs on the 
overall performance of the turbine. For more accurate 
computations, it may be necessary to include a local- 
ized 3-D analysis where flow reversals and flow stag- 
nation may cause the 1-D code predictions at these 
locations to be in error by a wide margin. 

NOMENCLATURE 
A area 

C p specific heat at constant pressure 

D diameter 

Dh hydraulic diameter defined by 

Du = 4A/P 

d pin fin diameter 

e surface roughness; or rib height 

f fanning friction factor defined by 

f - 6pD h /niV*p) 

Gr rotational Grashof number defined by 
Gr = N 7 r(l' w - T nw )D 3 0/v* 

II flow channel height 

h enthalpy 

h e coolant heat transfer coefficient 

j Colburn factor defined by 

j - h c rr'/*/( P VC„) 

k thermal conductivity 

t interval length 

m coolant mass flow rate at inlet to 

the channel 
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